Genetic regulation of volatile production in two melon introgression line collections with contrasting ripening behavior.

The importance of melon aroma in determining fruit quality has been highlighted in recent years. The fruit volatile profile is influenced by the type of fruit ripening. Non-climacteric fruits contain predominantly aldehydes, while climacteric fruits mainly produce esters. Several genes have been described to participate in volatile organic compounds (VOCs) biosynthesis pathways, but knowledge in this area is still incomplete. In this work we analysed the volatile profile of two reciprocal Introgression Line (IL) collections generated from a cross between 'Piel de Sapo' (PS) and 'Védrantais' (VED) melons, differing in their aroma profile and ripening behaviour. SPME GC-MS was performed to identify genes responsible for VOCs formation. More than 1000 QTLs for many volatiles were detected taken together both populations. Introgressions on chromosomes 3, 5, 6, 7 and 8 modified ester-aldehyde balance and were correlated to ripening changes in both genetic backgrounds. Some previously identified QTLs for fruit ripening might be involved in these phenotypes, such as ETHQV8.1 on chromosome 8 and ETHQV6.3 on chromosome 6. PS alleles on chromosomes 2, 6, 10 and 11 were found to increase ester content when introgressed in VED melons. Terpenes showed to be affected by several genomic regions not related to ripening. In addition, several candidate genes have been hypothesized to be responsible for some of the QTLs detected. The analysis of volatile compounds in two reciprocal IL collections has increased our understanding of the relationship between ripening and aroma and offers valuable plant material to improve food quality in melon breeding programs.

The haplotype-resolved T2T reference genome highlights structural variation underlying agronomic traits of melon.

Melon (Cucumis melo L.) is an important vegetable crop that has an extensive history of cultivation. However, the genome of wild and semi-wild melon types that can be used for the analysis of agronomic traits is not yet available. Here we report a chromosome-level T2T genome assembly for 821 (C. melo ssp. agrestis var. acidulus), a semi-wild melon with two haplotypes of ~373 Mb and ~364 Mb, respectively. Comparative genome analysis discovered a significant number of structural variants (SVs) between melo (C. melo ssp. melo) and agrestis (C. melo ssp. agrestis) genomes, including a copy number variation located in the ToLCNDV resistance locus on chromosome 11. Genome-wide association studies detected a significant signal associated with climacteric ripening and identified one candidate gene CM_ac12g14720.1 (CmABA2), encoding a cytoplasmic short chain dehydrogenase/reductase, which controls the biosynthesis of abscisic acid. This study provides valuable genetic resources for future research on melon breeding.

Regulation of climacteric fruit ripening in melon: recent advances and future challenges.

Fruit ripening is a complex and highly regulated process where tomato and strawberry have been the model species classically used for studying climacteric and non-climacteric fleshy fruit ripening types, respectively. Melon has emerged as an alternative ripening model because climacteric and non-climacteric cultivars exist, which makes it possible to dissect the regulation of ripening using a genetic approach. Several quantitative trait loci that regulate climacteric fruit ripening have been identified to date, and their combination in both climacteric and non-climacteric genetic backgrounds resulted in lines with different ripening behaviors, demonstrating that the climacteric intensity can be genetically modulated. This review discusses our current knowledge of the physiological changes observed during melon climacteric fruit ripening such as ethylene production, fruit abscission, chlorophyll degradation, firmness, and aroma, as well as their complex genetic control. From pioneer experiments in which ethylene biosynthesis was silenced, to the recent genetic edition of ripening regulators, current data suggest that the climacteric response is determined by the interaction of several loci under quantitative inheritance. The exploitation of the rich genetic diversity of melon will enable the discovery of additional genes involved in the regulation of the climacteric response, ultimately leading to breeding aromatic melon fruits with extended shelf life.

Thorough Characterization of ETHQB3.5, a QTL Involved in Melon Fruit Climacteric Behavior and Aroma Volatile Composition.

The effect of the QTL involved in climacteric ripening ETHQB3.5 on the fruit VOC composition was studied using a set of Near-Isogenic Lines (NILs) containing overlapping introgressions from the Korean accession PI 16375 on the chromosome 3 in the climacteric 'Piel de Sapo' (PS) genetic background. ETHQB3.5 was mapped in an interval of 1.24 Mb that contained a NAC transcription factor. NIL fruits also showed differences in VOC composition belonging to acetate esters, non-acetate esters, and sulfur-derived families. Cosegregation of VOC composition (23 out of 48 total QTLs were mapped) and climacteric ripening was observed, suggesting a pleiotropic effect of ETHQB3.5. On the other hand, other VOCs (mainly alkanes, aldehydes, and ketones) showed a pattern of variation independent of ETHQB3.5 effects, indicating the presence of other genes controlling non-climacteric ripening VOCs. Network correlation analysis and hierarchical clustering found groups of highly correlated compounds and confirmed the involvement of the climacteric differences in compound classes and VOC differences. The modification of melon VOCs may be achieved with or without interfering with its physiological behavior, but it is likely that high relative concentrations of some type of ethylene-dependent esters could be achieved in climacteric cultivars.

Fruit Morphology and Ripening-Related QTLs in a Newly Developed Introgression Line Collection of the Elite Varieties 'Védrantais' and 'Piel de Sapo'.

Melon is an economically important crop with widely diverse fruit morphology and ripening characteristics. Its diploid sequenced genome and multiple genomic tools make this species suitable to study the genetic architecture of fruit traits. With the development of this introgression line population of the elite varieties 'Piel de Sapo' and 'Védrantais', we present a powerful tool to study fruit morphology and ripening traits that can also facilitate characterization or pyramidation of QTLs in inodorous melon types. The population consists of 36 lines covering almost 98% of the melon genome, with an average of three introgressions per chromosome and segregating for multiple fruit traits: morphology, ripening and quality. High variability in fruit morphology was found within the population, with 24 QTLs affecting six different traits, confirming previously reported QTLs and two newly detected QTLs, FLQW5.1 and FWQW7.1. We detected 20 QTLs affecting fruit ripening traits, six of them reported for the first time, two affecting the timing of yellowing of the rind (EYELLQW1.1 and EYELLQW8.1) and four at the end of chromosome 8 affecting aroma, abscission and harvest date (EAROQW8.3, EALFQW8.3, ABSQW8.3 and HARQW8.3). We also confirmed the location of several QTLs, such as fruit-quality-related QTLs affecting rind and flesh appearance and flesh firmness.

Modulating climacteric intensity in melon through QTL stacking.

Fruit ripening is one of the main processes affecting fruit quality and shelf life. In melon there are both climacteric and non-climacteric genotypes, making it a suitable species to study fruit ripening. In the current study, in order to fine tune ripening, we have pyramided three climacteric QTLs in the non-climacteric genotype "Piel de Sapo": ETHQB3.5, ETHQV6.3 and ETHQV8.1. The results showed that the three QTLs interact epistatically, affecting ethylene production and ripening-related traits such as aroma profile. Each individual QTL has a specific role in the ethylene production profile. ETHQB3.5 accelerates the ethylene peak, ETHQV6.3 advances the ethylene production and ETHQV8.1 enhances the effect of the other two QTLs. Regarding aroma, the three QTLs independently activated the production of esters changing the aroma profile of the fruits, with no significant effects in fruit firmness, soluble solid content and fruit size. Understanding the interaction and the effect of different ripening QTLs offers a powerful knowledge for candidate gene identification as well as for melon breeding programs, where fruit ripening is one of the main objectives.

Knock-Out of CmNAC-NOR Affects Melon Climacteric Fruit Ripening.

Fruit ripening is an important process that affects fruit quality. A QTL in melon, ETHQV6.3, involved in climacteric ripening regulation, has been found to be encoded by CmNAC-NOR, a homologue of the tomato NOR gene. To further investigate CmNAC-NOR function, we obtained two CRISPR/Cas9-mediated mutants (nor-3 and nor-1) in the climacteric Védrantais background. nor-3, containing a 3-bp deletion altering the NAC domain A, resulted in ~8 days delay in ripening without affecting fruit quality. In contrast, the 1-bp deletion in nor-1 resulted in a fully disrupted NAC domain, which completely blocked climacteric ripening. The nor-1 fruits did not produce ethylene, no abscission layer was formed and there was no external color change. Additionally, volatile components were dramatically altered, seeds were not well developed and flesh firmness was also altered. There was a delay in fruit ripening with the nor-1 allele in heterozygosis of ~20 days. Our results provide new information regarding the function of CmNAC-NOR in melon fruit ripening, suggesting that it is a potential target for modulating shelf life in commercial climacteric melon varieties.

CRISPR/Cas9 gene editing uncovers the roles of CONSTITUTIVE TRIPLE RESPONSE 1 and REPRESSOR OF SILENCING 1 in melon fruit ripening and epigenetic regulation.

Melon (Cucumis melo) has emerged as an alternative model to tomato for studying fruit ripening due to the coexistence of climacteric and non-climacteric varieties. Previous characterization of a major quantitative trait locus (QTL), ETHQV8.1, that is able to trigger climacteric ripening in a non-climacteric background resulted in the identification of a negative regulator of ripening CTR1-like (MELO3C024518) and a putative DNA demethylase ROS1 (MELO3C024516) that is the orthologue of DML2, a DNA demethylase that regulates fruit ripening in tomato. To understand the role of these genes in climacteric ripening, in this study we generated homozygous CRISPR knockout mutants of CTR1-like and ROS1 in a climacteric genetic background. The climacteric behavior was altered in both loss-of-function mutants in two growing seasons with an earlier ethylene production profile being observed compared to the climacteric wild type, suggesting a role of both genes in climacteric ripening in melon. Single-cytosine methylome analyses of the ROS1-knockout mutant revealed changes in DNA methylation in the promoter regions of the key ripening genes such as ACS1, ETR1, and ACO1, and in transcription factors associated with ripening including NAC-NOR, RIN, and CNR, suggesting the importance of ROS1-mediated DNA demethylation for triggering fruit ripening in melon.

A cryptic variation in a member of the Ovate Family Proteins is underlying the melon fruit shape QTL fsqs8.1.

KEY MESSAGE: The gene underlying the melon fruit shape QTL fsqs8.1 is a member of the Ovate Family Proteins. Variation in fruit morphology is caused by changes in gene expression likely due to a cryptic structural variation in this locus. Melon cultivars have a wide range of fruit morphologies. Quantitative trait loci (QTL) have been identified underlying such diversity. This research focuses on the fruit shape QTL fsqs8.1, previously detected in a cross between the accession PI 124112 (CALC, producing elongated fruit) and the cultivar 'Piel de Sapo' (PS, producing oval fruit). The CALC fsqs8.1 allele induced round fruit shape, being responsible for the transgressive segregation for this trait observed in that population. In fact, the introgression line CALC8-1, carrying the fsqs8.1 locus from CALC into the PS genetic background, produced perfect round fruit. Following a map-based cloning approach, we found that the gene underlying fsqs8.1 is a member of the Ovate Family Proteins (OFP), CmOFP13, likely a homologue of AtOFP1 and SlOFP20 from Arabidopsis thaliana and tomato, respectively. The induction of the round shape was due to the higher expression of the CALC allele at the early ovary development stage. The fsqs8.1 locus showed an important structural variation, being CmOFP13 surrounded by two deletions in the CALC genome. The deletions are present at very low frequency in melon germplasm. Deletions and single nucleotide polymorphisms in the fsqs8.1 locus could not be not associated with variation in fruit shape among different melon accessions, what indicates that other genetic factors should be involved to induce the CALC fsqs8.1 allele effects. Therefore, fsqs8.1 is an example of a cryptic variation that alters gene expression, likely due to structural variation, resulting in phenotypic changes in melon fruit morphology.

A novel introgression line collection to unravel the genetics of climacteric ripening and fruit quality in melon.

Introgression lines are valuable germplasm for scientists and breeders, since they ease genetic studies such as QTL interactions and positional cloning as well as the introduction of favorable alleles into elite varieties. We developed a novel introgression line collection in melon using two commercial European varieties with different ripening behavior, the climacteric cantalupensis 'Védrantais' as recurrent parent and the non-climacteric inodorus 'Piel de Sapo' as donor parent. The collection contains 34 introgression lines, covering 99% of the donor genome. The mean introgression size is 18.16 Mb and ~ 3 lines were obtained per chromosome, on average. The high segregation of these lines for multiple fruit quality traits allowed us to identify 27 QTLs that modified sugar content, altered fruit morphology or were involved in climacteric ripening. In addition, we confirmed the genomic location of five major genes previously described, which control mainly fruit appearance, such as mottled rind and external color. Most of the QTLs had been reported before in other populations sharing parental lines, while three QTLs (EAROQP11.3, ECDQP11.2 and FIRQP4.1) were newly detected in our work. These introgression lines would be useful to perform additional genetic studies, as fine mapping and gene pyramiding, especially for important complex traits such as fruit weight and climacteric ripening.

Genetic dissection of aroma biosynthesis in melon and its relationship with climacteric ripening.

Aroma is an essential trait in melon fruit quality, but its complexity and genetic basis are still poorly understood. The aim of this study was the identification of quantitative trait loci (QTLs) underlying volatile organic compounds (VOCs) biosynthesis in melon rind and flesh, using a Recombinant Inbred Line (RIL) population from the cross 'Piel de Sapo' (PS) × 'Védrantais' (VED), two commercial varieties segregating for ripening behavior. A total of 82 VOCs were detected by gas chromatography-mass spectrometry (GC-MS), and 166 QTLs were identified. The main QTL cluster was on chromosome 8, collocating with the previously described ripening-related QTL ETHQV8.1, with an important role in VOCs biosynthesis. QTL clusters involved in esters, lipid-derived volatiles and apocarotenoids were also identified, and candidate genes have been proposed for ethyl 3-(methylthio)propanoate and benzaldehyde biosynthesis. Our results provide genetic insights for deciphering fruit aroma in melon and offer new tools for flavor breeding.

QTLs and candidate genes analyses for fruit size under domestication and differentiation in melon (Cucumis melo L.) based on high resolution maps.

BACKGROUND: Melon is a very important horticultural crop produced worldwide with high phenotypic diversity. Fruit size is among the most important domestication and differentiation traits in melon. The molecular mechanisms of fruit size in melon are largely unknown. RESULTS: Two high-density genetic maps were constructed by whole-genome resequencing with two F2 segregating populations (WAP and MAP) derived from two crosses (cultivated agrestis × wild agrestis and cultivated melo × cultivated agrestis). We obtained 1,871,671 and 1,976,589 high quality SNPs that show differences between parents in WAP and MAP. A total of 5138 and 5839 recombination events generated 954 bins in WAP and 1027 bins in MAP with the average size of 321.3 Kb and 301.4 Kb respectively. All bins were mapped onto 12 linkage groups in WAP and MAP. The total lengths of two linkage maps were 904.4 cM (WAP) and 874.5 cM (MAP), covering 86.6% and 87.4% of the melon genome. Two loci for fruit size were identified on chromosome 11 in WAP and chromosome 5 in MAP, respectively. An auxin response factor and a YABBY transcription factor were inferred to be the candidate genes for both loci. CONCLUSION: The high-resolution genetic maps and QTLs analyses for fruit size described here will provide a better understanding the genetic basis of domestication and differentiation, and provide a valuable tool for map-based cloning and molecular marker assisted breeding.

Genetic dissection of climacteric fruit ripening in a melon population segregating for ripening behavior.

Melon is as an alternative model to understand fruit ripening due to the coexistence of climacteric and non-climacteric varieties within the same species, allowing the study of the processes that regulate this complex trait with genetic approaches. We phenotyped a population of recombinant inbred lines (RILs), obtained by crossing a climacteric (Védrantais, cantalupensis type) and a non-climcteric variety (Piel de Sapo T111, inodorus type), for traits related to climacteric maturation and ethylene production. Individuals in the RIL population exhibited various combinations of phenotypes that differed in the amount of ethylene produced, the early onset of ethylene production, and other phenotypes associated with ripening. We characterized a major QTL on chromosome 8, ETHQV8.1, which is sufficient to activate climacteric ripening, and other minor QTLs that may modulate the climacteric response. The ETHQV8.1 allele was validated by using two reciprocal introgression line populations generated by crossing Védrantais and Piel de Sapo and analyzing the ETHQV8.1 region in each of the genetic backgrounds. A Genome-wide association study (GWAS) using 211 accessions of the ssp. melo further identified two regions on chromosome 8 associated with the production of aromas, one of these regions overlapping with the 154.1 kb interval containing ETHQV8.1. The ETHQV8.1 region contains several candidate genes that may be related to fruit ripening. This work sheds light into the regulation mechanisms of a complex trait such as fruit ripening.

Transposons played a major role in the diversification between the closely related almond and peach genomes: results from the almond genome sequence.

We sequenced the genome of the highly heterozygous almond Prunus dulcis cv. Texas combining short- and long-read sequencing. We obtained a genome assembly totaling 227.6 Mb of the estimated almond genome size of 238 Mb, of which 91% is anchored to eight pseudomolecules corresponding to its haploid chromosome complement, and annotated 27 969 protein-coding genes and 6747 non-coding transcripts. By phylogenomic comparison with the genomes of 16 additional close and distant species we estimated that almond and peach (Prunus persica) diverged around 5.88 million years ago. These two genomes are highly syntenic and show a high degree of sequence conservation (20 nucleotide substitutions per kb). However, they also exhibit a high number of presence/absence variants, many attributable to the movement of transposable elements (TEs). Transposable elements have generated an important number of presence/absence variants between almond and peach, and we show that the recent history of TE movement seems markedly different between them. Transposable elements may also be at the origin of important phenotypic differences between both species, and in particular for the sweet kernel phenotype, a key agronomic and domestication character for almond. Here we show that in sweet almond cultivars, highly methylated TE insertions surround a gene involved in the biosynthesis of amygdalin, whose reduced expression has been correlated with the sweet almond phenotype. Altogether, our results suggest a key role of TEs in the recent history and diversification of almond and its close relative peach.

Mapping Cucumber Vein Yellowing Virus Resistance in Cucumber (Cucumis sativus L.) by Using BSA-seq Analysis.

Cucumber vein yellowing virus (CVYV) causes severe yield losses in cucurbit crops across Mediterranean countries. The control of this virus is based on cultural practices to prevent the presence of its vector (Bemisia tabaci) and breeding for natural resistance, which requires the identification of the loci involved and the development of molecular markers for linkage analysis. In this work, we mapped a monogenic locus for resistance to CVYV in cucumber by using a Bulked Segregant Analysis (BSA) strategy coupled with whole-genome resequencing. We phenotyped 135 F3 families from a segregating population between a susceptible pickling cucumber and a resistant Long Dutch type cucumber for CVYV resistance. Phenotypic analysis determined the monogenic and incomplete dominance inheritance of the resistance. We named the locus CsCvy-1. For mapping this locus, 15 resistant and 15 susceptible homozygous F2 individuals were selected for whole genome resequencing. By using a customized bioinformatics pipeline, we identified a unique region in chromosome 5 associated to resistance to CVYV, explaining more than 80% of the variability. The resequencing data provided us with additional SNP markers to decrease the interval of CsCvy-1 to 625 kb, containing 24 annotated genes. Markers flanking CsCvy-1 in a 5.3 cM interval were developed for marker-assisted selection (MAS) in breeding programs and will be useful for the identification of the target gene in future studies.

A comprehensive genome variation map of melon identifies multiple domestication events and loci influencing agronomic traits.

Melon is an economically important fruit crop that has been cultivated for thousands of years; however, the genetic basis and history of its domestication still remain largely unknown. Here we report a comprehensive map of the genomic variation in melon derived from the resequencing of 1,175 accessions, which represent the global diversity of the species. Our results suggest that three independent domestication events occurred in melon, two in India and one in Africa. We detected two independent sets of domestication sweeps, resulting in diverse characteristics of the two subspecies melo and agrestis during melon breeding. Genome-wide association studies for 16 agronomic traits identified 208 loci significantly associated with fruit mass, quality and morphological characters. This study sheds light on the domestication history of melon and provides a valuable resource for genomics-assisted breeding of this important crop.

An Improved Melon Reference Genome With Single-Molecule Sequencing Uncovers a Recent Burst of Transposable Elements With Potential Impact on Genes.

The published melon (Cucumis melo L.) reference genome assembly (v3.6.1) has still 41.6 Mb (Megabases) of sequences unassigned to pseudo-chromosomes and about 57 Mb of gaps. Although different approaches have been undertaken to improve the melon genome assembly in recent years, the high percentage of repeats (~40%) and limitations due to read length have made it difficult to resolve gaps and scaffold's misassignments to pseudomolecules, especially in the heterochromatic regions. Taking advantage of the PacBio single- molecule real-time (SMRT) sequencing technology, an improvement of the melon genome was achieved. About 90% of the gaps were filled and the unassigned sequences were drastically reduced. A lift-over of the latest annotation v4.0 allowed to re-collocate protein-coding genes belonging to the unassigned sequences to the pseudomolecules. A direct proof of the improvement reached in the new melon assembly was highlighted looking at the improved annotation of the transposable element fraction. By screening the new assembly, we discovered many young (inserted less than 2Mya), polymorphic LTR-retrotransposons that were not captured in the previous reference genome. These elements sit mostly in the pericentromeric regions, but some of them are inserted in the upstream region of genes suggesting that they can have regulatory potential. This improved reference genome will provide an invaluable tool for identifying new gene or transposon variants associated with important phenotypes.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops.

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a 'SyntenyViewer' to view genome synteny between different cucurbit species and an 'RNA-Seq' module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.